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ABSTRACT 

Two computer simulation models of earthquakes are studied for the depend- 
ence of the pattern of events on the model assumptions and input parameters. 

Both models are adaptations from the work of Dieterieh (Dieterich, 1972), and 
represent the seismically active region by mechanical blocks which are connected 
to one another and to a driving plate. The blocks slide on a friction surface. In 
the first model we employ elastic forces and time independent friction to simulate 
main shock events. We find that the size, length, and time and place of event 
occurrence are influenced strongly by the magnitude and degree of homogeniety in 
the elastic and friction parameters of the fault region. For example, periodically 
reoccurring similar events are frequently observed in simulations with near 
homogeneous parameters along the fault, whereas, seismic gaps are a common 
feature of simulations employing large variations in the fault parameters. The 
second model incorporates viscoelastic forces and time-dependent friction to 
account for aftershock sequences. The periods between aftershock events in- 
crease with time and the aftershock region is confined to that which moved in 
the main event. 
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I. INTRODUCTION 

In the absence of a general theory of motion along a fault there are at least 
five avenues of approach to understanding earthquake dynamics. These avenues 
Include field Investigations, laboratory experiments, simplified theoretical an- 
alyses, model building, and computer-numeric simulations. The present paper 
focuses on a computer simulation. There are several reasons why such simula- 
tions are attractive research tools. The first of these is time compression. 
Events spanning many decades, even centuries, can be simulated in a few sec- 
onds on a high speed computer. The second is convenience in investigating al- 
ternative hypotheses. All the conditions of the simulation are under the control 
of the progr mmer and the differing consequences of alternative assumptions are 
readily identified. The third attractive feature is the great detail with which the 
motton can be observed in the simulation. This can be a great aid in explaining 
the consequences of the motion. The major difficulty with computer simulations 
is the question of correspondence with "real world" phenomena. Highly idealized 
models, which are clear in conceptual detail, may not be adequate representations 
of naturally occurring situations. On the other hand complex simulations are lim- 
ited by our lack of knowledge about naturally occurring conditions and by restric- 
tions on computer speed and storage requirements. Furthermore, complex 
simulations may be so confusing in detail that they are uninstructive in basic 
causes and effects. In the present effort we are aiming primarily for physical 
clarity and nave, therefore, used conceptually simple ideas and techniques. 
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Despite this limitation, much of the behavior exhibited by our simulator mimics 
natural behavior and certain observed natural phenomena find ready explanation 
through this technique. 

Most simulators are based on the model of coupled massive blocks first in- 
troduced by Burridge and Knopoff (Burridge and Knopoff, 19(57) and illustrated 
in Figure 1. In this model 'he fault is divided into coupled blocks which slide on 
a frictional surface. The coupling which corresponds to the elastic and viscous 
properties of the material is represented by various combinations of springs and 
dashpots. The blocks are driven through a coupling to a moving plate. The plate 
moves with a drift velocity, u, which is unaffected by the frictional resistance along 
the fault. The one-dimensional model of Burridge and Knopoff included ten blocks 
some with predominate elastic coupling lo their neighbors and some with pre- 
dominate Maxwell-viscous coupling. Dletorich (Dieterich, 1972) elaborated on 
the model by including fifty blocks and using the standard linear solid coupling 
shown in Figure 1. Otsuka (Otsuka, ’972) developed a two-dimensional simula- 
tor with elastic couplding, and Dieterich (Dieterich, 1973) reported on a three- 
dimensional, elastic coupling, simulator. 

In the present work we will examine in more detail some one-dimensional 
simulators, primarily those due to Dieterich. We will examine, in particular, 
such questions as how the pattern of stress release in simulator earthquakes 
(SEQe) depends on the frictional characteristics of the fault and how simulator 
aftershock sequences depend on material properties. 
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II. ELASTIC MODEL 


The simplest model that we have employed assumes perfect elastic coupling 
between a block and its nearest neighbors and to the driving plate. Consider the 
elastic forces acting on the ith block of Figure 1 

Fj * Kj(X| + , - Xj) + . | (Xj . , - Xj ) + Kj'iut - Xj) (1) 

where u is the velocity of the driving plate, t is the time, Xj * , , Xj _ ,, and Xj 
are the displacements of the 1 + 1, i - 1, and ith blocks, K, and K j _ , are the 
spring constants of the i and i - 1 th connecting springs and K’ is the spring con- 
stant of the leaf spring connecting block i to the driving plate. To each block we 
also associate a static frictional strength, f ] , which resists motion. Once the 
pulling force overcomes the frictional resistance to motion the dynamics of the 
ith block is governed by the force equation 


ntjXi = Fi - fj 

= Kj(Xj + | - Xj) + Kj _ , (Xj _ , - x ( ) + Kj'(ut — Xj> — fj* 


(2) 


where nij is the block mass and fj is the dynamic friction force. 


Consider now the situation with all blocks initially at rest and F, < f^ for 
each block. As time advances F ( increases due to the motion of the driving block 
and the consequential increase in the driving spring tension. We can associate 
with each block a time, t", at which time the block woi Id first move if unaffected 
by the motion of any other block. Thus 


t? = - 


K; K. , „ if 

— (Xj 1 - , ) + (X? - X? . ) + — 

I/I 1 1*1 I/I I II I/I 


( 3 ) 
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That block which has the minimum value of tj' Is the first to move. The subse- 
quent motion of all blocks can be determined as follows: No block moves until 
Fj f* for that block. At that time the block's initial acceleration Is 

fj> - n* if - rf 

••<> _ t i _ i i . . 


The subsequent dynamics is dictated by Equation 2. Initially the block's acceler- 
ation Is positive and Its position and velocity increase. As the tension on the 
block is relieved >y the sliding motion the acceleration decreases Incoming neg- 
ative for sufficient spring compression. The velocity begins to decrease and 
ultimately goes to zero :*t which point motion of that block terminates if IF. I 
< f*. It is evident from proceeding equations, however, that the motion of 
one clock will cause compression or expansion of the connecting springs to its 
neighbors. This can stimulate motion of the adjoining blocks thereby propagating 
the earthquake along the fault. Thu.? the motion along the fault is governed by a 
set of coupled dynamic equations of the form of Equation 2. Numeric techniques 
for solving such equations are well known ami we need comment here only that 
the solution would be ambiguous if we did not impose some boundary condition of 
the end springs. In the simulations reported here we have assumed periodic 
boundary conditions for N blocks (N = 50 in most of our simuitations). 


Before continuing the discussion it is instructive to solve analytically the 


simple case in which only one block moves during an earthquake. We further 
assume K ( = K t + , = Kj and take as initial conditions x ■' - xj' + , x" = 0. Then 






suppressing the subscript 1 


.. 3K Kut-P 

X ♦ — x * 

m m 


Since the velocity of the driving plate Is very slow (~5cm/yr) the quantity Kut 
is essentially unchanged during the motion and we let t *= t„ constant. "njus the 


equations of motion becomes 


.. 3K Kut 0 - P 
x + — x * 

m m 


which has the simple solution 


2 / f d \ , , /Jk (t-i 

7 ( u, °.-)s,n- /- — 


But the condition for starting motion is Kut 0 • f s so that 

2 f* , /IK 1 f d 

x a - — (1 - a)sin z / — - ; a = — 

3 k V m 2 f* 


We take note of several easily derived consequences. First the tots' displace- 
ment of block i due to this event is 

2 f* 

*m.x « 


while the duration of the event is 


•nux " *o 


Notice that the duration of the single block event depends only or. the mass of the 


block and the spring constant and is unaffected by the friction parameters. The 
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peak velocity is 


(f»-f d ) 

\/TKm 




The typical time scale for re-occurrence of motion in this block is 


AT - 


f*-f d 


The Initial force on block 1 was Kut„. The final force is ^-K(ut 0 - 3x max ). Thus 
the change in the stress force, AF, is 

AF i " 3K *max * fd) < 13 > 


2(1 - a) 


The stress force rise in the adjacent block (e.g. , AF, + , ) is 

AF i* i * * \ (f, - fd > < 15 > 

AF ,*i o 

-ar--r(i-a) < 16 > 

F l ♦ l 3 

Notice that the fractional stress changes depend only on the ratio of dynamic to 
static friction, independent of other material properties. If, for example, a 
0.8, there is a forty percentdrop in the stress of the ith block due to this event and 
a thirteen percent increase in the stress of its neighbors. If this additional stress 
rise were sufficient to overcome the frictional strength of the neighbors these 
blocks would have been stimulated into motion. Furthermore, the absolute change 
in the stress force as a result of the earthquake is determined by twice the 
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difference between the static and dynamic frictional strengths, again independ- 
ent of other material parameters. 

III. VISCOELASTIC MODEL WITH TIME-DEPENDENT FRICTION 

The clastic model which we discussed In the last section generates simula- 
tions which have many features common to naturally occurring events. These 
features include the slow accumulation of tectonic .stress and its rapid release by 
stick-slip sliding in earthquakes, the existence of seismic gaps in areas on high 
frictional and elastic constant variability, and the gradual propagation of a se- 
cession of earthquake events through regions of relative uniform properties. 

I 

Events of varying magnitude, ground displacements, and volume are generated. 
Nevertheless there arc a number of important phenomena which cannot l»e simu- 
lated with this model. Foremost among these are aftershock sequences. In or- 
der 10 generate aftershocks we assume some mechanism for the rapid recovery 
of all or part of the stress released by the main shock is necessary. The tectonic 
forces represented by the slowly stretching clastic springs build too slowly to 
provide an aftershock mechanism. By contrast, Burrldge and Knopoli (Bur ridge 
and Knopoff, 1967) showed that viscous creep following an earthquake can redis- 
tribute stress in a manner which produces aftershocks. Dieterich (Dieterich, 
1972) proposed partial stress recovery due to viscoelasticity in a region which 
slides during an earthquake then remains locked while the viscous forces adjust 
to the displacement. When coupled with a post-earthquake weakening oi the fault 




frictional strength this mechanism will also generate aftershock sequences. We 


now elaborate on this model. 


Consider the same mechanical blocks as in the previous section. Now, how- 
ever, instead of having coupling by simple elastic springs, the coupling is a par- 
allel combination of elastic and Maxwellian elements as shown in insert b of Fig- 
ure 1. Assuming thut the blocks are not moving and that the driving plate velo- 
city is still u. The viscoelastic force equation is 


♦ K^X,* , -XjHI -c" 0 "‘•♦ )/T ») + K, . ,(x, . , -x,)(l - e ' ( ‘ - ') 


+ K', urj(l -c -< * ♦ Kju(t - tjC~^ 1 " ,,)/T| ) 

- K!x,(l -c 0 ■'• )/r ‘) 


(17) 


'0 


*!♦ 


In the proceeding equation we have introduced the following symbols: 

Gj - force on spring i at time t^ 

Gj . ( = force on spring i - 1 at time tj_ 

force on spring i’ at time t, 
last time block i or block i + 1 moved 

last time block i or block 1 - 1 moved 

t, - last time block i moved. 

The tau's (r) are the relaxation times of the viscoelastic system. The condition 
for the onset on a simulation earthquake is 

Fj(t) * fj(t) (18) 


t 


i- 
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where the atulic frictional strength f*(t), is also a weak function of time increas- 
ing slowly from some minimum value immediately following cessation of block 
motion. Equation 17 Is a transcendental equation for which can be solved by anv 
of several numeric techniques. The acceleration and de-acceleration occurring 
during an earthquake occur on a time scale much faster than the viscous response 
time hence we can Ignore viscous properties during the earthquake and write the 
following equation of motion 

m'Xj = [K t (x, , - Xj)+ Kj _ , (x t _ , - Xj) ♦ K'tut - x,)] 

+ [K : i(x, + , -x,)-(x® + , -xf)) + K 2 ((Xj _ , - Xj) - (xp_ , - x”)) (19) 

+ k; Out - Xj ) - (utjj - x?))'] + |3” + 3? + +3“’ I -fj d 
In this equation the first bracketed quantity represents the elastic force due to 
the stretching of the purely elastic element of Figure lc. The second and third 
bracketed quantities represent the forces on the Maxwell elements. The second 
bracket gives ne chang in the elastic force due to stretching or compression of 
Maxwell elements springs from their position at beginning of simulation earth- 
quake. The third bracket gives the force prior to this stretching or compression. 
This model has two key elements: the viscoelasticity is responsible for a partial 
recovery of the stress drop occurring during an earthquake event and the time 
dependent friction is responsible for a weakening of the fault strength. The model 
operates as follows. Assume initially (for simplicity only) that all elastic and 
Maxwell elements are relaxed, that there is no stress acting on any of the blocks. 
As in the purely elarMo mode 1 , stress accumulates with the drift of the driving 



plate until one block Is stressed beyond Its frictional strength and begins to 
slide. The subsequent motion may involve one or more blocks depending on the 
frictional strength of the other blocks. The motion of the blocks results in a 
stress drop whose magnitude is determined by the sum of the products of the 
spring constants and the change in spring stretch. This is all similar to the pure 
elastic case. Now, however, while the blocks remain fixed following the simu- 
lated earthquake the viscous dashpots tend to aejust to relax the stretch or com- 
pression in the Maxwell element springs. Consider the primed coupling element 
which has undergone a reduction in the spring stretch AL. The resultant stress 
drop is (K' + K',) AL. If the Maxwell element dashpot has sufficient time to re- 
lax the change in the stretch of the K', spring, then the final stress drop following 
viscous adjustment is K' AL. Thus the fractional stress reco\ery of the original 
stress drop is K,/(K' + K' : ). Assume, lor illustrative purposes that Kj = 3K' , 
then seventy five percent of the original stress drop durii.g a simulation earth- 
quake is recovered during the subsequent stress rise due to viscous adjustment. 

If the frictional strength of the fault in this region is sufficiently reduced after the 
ear*h 7 »ake compared to the pre-earthquake level and if the viscous response is 
sufficiently rapid, then an aftershock will be generated. The process may be re- 
peated many times. The reduction in frictional strength along the displaced por- 
tion of the fault is a consequence of the finite time required to heal the break sub- 
sequent to the slippage. This model uses a friction strength which rises loga- 
rithmically with time sLnce the last slip. The time scai° for the aftershocks is 


10 


T 


I 


T 


usually of the* same order as the viscous relaxation time. We will present il- 
lustrations of these aftershock sequences in the following section. We should 
emphasize that this aftershock theory is just one of several alternative treat- 
ments of the aftershock problem. Subsequent papers will deal with computer 
simulations employing other theoretical models. 

IV. RESULTS 

W r e now discuss the results of simulations using the models discussed in the 
previous sections. We first consider the simple case of purely elastic coupling 
between the blocks and the driving plate. Figure 2a shows block acceleration, 
velocity, and displacement versus time for typical small two block event. This 
SEQ was initiated when the elastic forces acting on block 44 overcame the static 
frictional strength of the block. The block begins to move with a pea* accelera- 
tion 11 em/sec : . * As this block slides to the right and begins to relax the force 
acting on it, additional force is applied to blocks 43 and 45 by the compression 
and extension of the i - 43 and i - 44 springe. At ~7 seconds after the initiation 
of the event, the force on block 43 has reached the frictional strength limit of that 
bio and it too begins to move to the right. As both blocks continue to slide to 
the right additional force is applied to blocks 42 and 45. Nevertheless these 
blocKS remain unmoved as even the additional stress due to the motion of blocks 

*The acceleration, velocity, displacement, and time axis shown on Figure 2 are scaled by the independent 
parameters f, K, and m. Thus, for example, the forty second time scale of Figure 2 is arbitrary to die 
extent that independent parameters can be changed. Sec the Appendix for the scale relationships. 
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•14 und 43 is insufficient to overcome the frictional strengths. The moving blocks 
continue to de accelerate until they come to rest at ~34 und 39 seconds after the 
initiation of the event. The force conditions on the fault proceeding and following 
this SEQ are also shown in Figure 2. Notice the stress relaxation in the moved 
blocks and the stress buildup in the adjacent unmoved blocks, in the case where 
the frictional and elastic properties of the fault region vary only slightly from 
block to block, this predisposes this unmoved border region to motion in a sub- 
sequent SEQ. By this mechanism the location of SEQ'b propagates down the en- 
tire fault leaving few, if any, seismic gaps. We can see this more clearly in an 
event plot as shown in Figure 3. In this plot each line summarizes pertinent in- 
formation concerning one SEQ. The first column gives the event number, the 
second the time of occurrence of the SEQ. The third column of fifty spaces rep- 
resents the fifty blocks of the fault used in our simulations. Motion of a particu- 
lar block during the event is indicated by an X; otherwise the space is left blank. 
Thus in event 316 at time ^5. 147 x 1 0 1 " seconds, blocks 6 through 9 moved and 
all other blocks remained stationary. The fourth column gives the block number 
of the first block to move in the event, and the ftfth column tells the tota* number 
of blocks displaced in the SEQ. Figure 3 is an event plot for which all K’ = 1 
x 10 1 dyne/cm, Kj = Kj/2, and the f'J were in the range (1.95-2.05)x 10 2 ° dynes. 
All the blocks have an equal mass of 2.8 x 10 2h grams. Notice that in this case 
there is a cycle of similar events which tend to repeat themselves. Although 
there is no true periodicity of events, there Is an approximate characteristic 
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time interval of about 2. 5 x 10 h seconds during which almost all blocks along the 
fault move at least once. 

By contrast to this case where the elastic and frictional properties vary only 
slightly from block to block, wc now consider a case. Figure 4, in which the 
friction is allowed to varv randomly from 1 x 10*’" to 3 x 10 : " dynes. In this 
case there is considerably more scatter in the location of SKQ's. 

Figures 5 and 6 show two other cases of some interest. In Figure 5 the 
frictional strength varies smoothly from a minimum value of about 1 x 10 ?0 
dynes at blocks 1 and 50 to a maximum value near 3 x 10' 1,1 dvnes in the middle 
of the fault. The function chosen to represent this variation is f^ 3. 16 x 10*’" 
expMi - 25. 5)/24. 5) " - 0. 16 x 10*’" (dvnes). Because of the symmetry of the 
fault equal simultaneous events occur on the right and left sides of the figure. 

Two distinct patterns are observed. In the early stages of the simulation the 
pattern of SEQ events traces out the friction curve. The weaker elements move 
frequently with small slippages while the stronger elements are more resistant 
to being displaced. Gradually, however, the pattern evolves into one in which 
large sections of the fault move in discrete steps separated by intervals of little 
activity in those sections. The reoccurrence time between the large events de- 
creases with time while the number of blocks involved increases. After the larg- 
est of these events has been completed much o f the tension in these sections of 
the fault has been relieved, and the earlier pattern of tracing out the friction 
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curve is repeated. This behavior is characteristic of a situation in which the 
fault goes through a charge -discharge cycle. During the charging phase stress 
gradually accumulates and the small events which do occur onl> partially relieve 
the mounting stress. Anally a large event or events discharge most of the ac- 
cumulated stress after which the cycle begins anew. 

In Figure 6 the friction is, as in Figures 3 and 4, randomly distributed, but 
now the limits on the random number generator are 0-1 x 10 : " dynes. The par- 
ticular distribution used is shown on the figure. We notice in Figure fi that the 
weak elements move frequently, often in single block events. 

A critical factor in determining the number of blocks that move in a given 
SEQ is the relative magnitudes of the connecting and leaf springs. Larger values 
of the leaf spring constants (relative to that of the connecting springs) favor events 
with a large number of blocks moving. The explanation is simply that with the 
larger leaf spring constant the tension between the moving plate and the associated 
connected block can be relieved with a relatively small displacement. This small 
displacement coupled to the weaker spring constant of the connecting spring is 
less likely to induce motion in the adjacent block than it would if the situation 
were reversed. 

Three examples of simulations using the viscoelastic model arc shown in 
Figure 7. The effect of the longer time constant in reducing the number of after- 
shocks and in increasing their inter-event period is obvious. Results generated 
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by the viscoelastic model are very sensitive to the relative values of the elastic 


and viscous parameters. We have performed simulations with many combinations 

K' 

of these parameters in the ranges 10 16 dynes <K' < 10 dynes, 0. 1 < 10, 


k; 




1. 5 < — and — < 5, 10 < r\ and t «S 10 4 , 0 < f* < 3 x 10 20 dynes and 0 < — 
K K' f» 


< 0. 99. 


Several general observations can be made. 

1. For simulations in which K K' and K 2 K 2 the number of blocks mov- 
ing during a typical event rises as the system advances from an unstressed to 
highly stressed state. Eventually a condition is reached when a high fraction of 
the simulation events involve all or nearly all blocks. This condition is influenced 
to some extent by the choice of periodic boundary conditions. 

2. Conditions favoring a large number of aftershocks include small visco- 
elastic response times and larger values of K' 2 with respect to K'. The former 
condition assures that the full effects of the viscoelastic stress recovery will be 
rapid compared to the time for appreciable rise in the frictional strength. The 
second condition implies that an appreciable portion of the initial stress drop 
during the SEQ is subsequently reduced so that the post-SEQ stress level ap- 
proaches the pre-SEQ level. 

3. The aftershocks tend to have smaller displacements than the main 
shock. This is consistent with the trend observed in nature and with the fact that 
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only a partial stress recovery occurs. It is not true however, that there Is a 
monotone decrease in the average displacement with aftershock number. The 
exact nature of the aftershock displacements depends in a complicated manner on 
the frictional characteristics, viscoelastic parameters, and number of blocks in- 
volved in the event. 

4. The mean time between aftershocks tends to increase with time. For 
som ? simulations the time interval between the main and first aftershock is the 
mitimum time interval between events in a sequence. In other cases the mini- 
mum time interval occurs between the first and second aftershocks. In either 
case there is a subsequent trend tow'ard increasing periods between aftershocks 
provided aftershocks with a common epicentral block are considered. For 
aftershocks with different epicentral blocks there may be a super-position of 
several sequences, each involving its own epicentral block. Alternatively, the 
behavior may be more complicated. Some of the simulated aftershock sequences 
show a logarithmic increase in the number of aftershocks with time, consistent 
with Omori's (Omori, 1894) observation of a (a + bt) _l dependence of aftershock 
frequency on time following the main shock. In other cases there are consider- 
able deviations from this simple tune dependence. 

5. The aftershock region is generally confined to the region of the main 
shock. The number of blocks displaced in an aftershock is less than or some- 
times equal to the number displaced in the main event. 
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V. CONCLUSIONS 

In this paper we have begun our examination ol various numerical simulators 
of earthquakes. We have focused attention on two theoretically simple models, 
one employing elastic forces and time -independent friction, the other employing 
viscoelastic forces and time -dependent friction. Our central conclusions are as 
follows : 

1. When the clastic and frictional parameters are nearly uniform from 
block to block the SKQ's are characterized by a propagation of successive events 
down the fault. Conversely, seismic gaps ure more likely to occur when there 
are large variations in the elastic or frictional parameters. 

2. The periodic reoccurrence of similar SEQ events is frequently observed 
in cases where the parameters of the fault and the blocks are nearly homogeneous 
from location to location. The reoccurrence time is determined primarily by the 
difference between the static and dynamic friction and by the elastic constants of 
the springs. 

3. Simulator events involving large displacements tend to involve the mo- 
tion of several blocks. The simultaneous motion of adjacent blocks lessens the 
displacement restricting elastic forces between blocks. 

4. Frictionally weak elements tend to move frequently with cone- mitant 
small average displacements. 
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K* 

5. As the ratio, — , of the elastic driving force constant to connecting 

K 

force constant increases, the length of the average SfcXi event decreases. 

6. In the viscoelastic model the aftershocks exhibit several characteristic 
features: 

a. the period between aftershocks increases with time following the 
main event, 

b. the magnitudes and block displacements of the aftershocks are less 
than that of the main shock. However, there is no consistent de- 
crease in aftershock magnitude with time, 

c. the aftershocks are confined to the region of the main event. Fewer, 
or at most the same, number of blocks are involved in aftershocks. 

7. The occurrence of aftershocks in the viscoelastic model is enhanced by 

n 

short viscoelastic time constant, r, and by large values of the ratio 

Several features concerning numerical simulators await further investiga- 
tion. These features include, among others, the exploration of alternative earth- 
quake theories particularly with respect to fault friction instability and aftershock 
occurrence. Also deserving examination in the correlation among the inj>ut and 
output parameters of the model. For example. King (King, 1975) has observed 
in his mechanical models the statistically valid relationship: log (Energy released) 
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* 


■ A + B log(Nd) where N Is the number of blocks Involved In the event, <1 is the 
average block displacement, and B*l, We have observed a similar correlation 
in our numerical simulations. Finally a detailed examination of the im|K>rtance 
of boundary conditions and geometry and topology in two and three dimensional 
models also awaits further discussion. We hope to address some of thrBc ques- 
tions in subsequent papers. 
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FIGURE CAPTIONS 


Figure 1. Mechanical block representation of fault. Each block which rests on 
the friction surface is coupled to nearest neighbors and to the driving 
plate which moves to the right with velocity u. (a) elastic coupling ele- 
ment - perfect spring with spring constant K. (b> viscoelastic coupling 
with spring constants K and K 2 and dashpol viscosity 17 . 

Figure 2. Example of a two block SEQ. The displacement velocity, ami acceler- 
ation of blocks 13 and 41 are shown as a function of time. This SEti 
occurred early in a simulation when the pre-SEl^ stresses on blocks 
42-45 were equal. After the SEt^ the stress on the moved blocks has 
been reduced while the stress on the adjacent unmoved blocks has 
risen. 

Figure 3. Event plot for simple elastic force model with K[ 1 x 10 1 dyne/cm, 

K, = 5 x lO 1 ' 1 dyne/cm and f^ randomly distributed in r. .igc (2 t 0 . 5) x 
10 20 dyne. 

Figure 4. Event plot for simple elastic force model with Kj and K, as in Figure 
3 but fj 1 randomly distributed in range (2 1 1) x 10 :o dyne. 

Figure 5. Event plot for simple elastic force model with K,’ and K, as in Figures 
3 and 4 and f* « 3. 1(5 x 10 20 exp(-(i - 25. 5)/24. 5 ) 2 - 0. 16 x 10 20 


(dynes). 


FIGURE CAPTIONS (Continued) 


Figure 6. Event plot for simple elastic force model with K[ and K, as in Figures 
3-5 but f* randomly distributed in range 0 < f' < 1 x 1 0 ’ ' dyne. 

Figure 7. Event plot for viscoelastic force model with time dependent friction. 
The parameters of the simulation are K{ - 1 x 10 U , Kj 3K, f K, 

- 1 x 10' 7 , K 2( - 3K,. lx 10 20 dyne < f‘(t 0 ) < 3 x 10 20 dyne, (■ 
-0.98f‘(t 0 ), f[(t)-fj(t 0 )(l * 0.021 log (1 + t-t 0 )|. 



40000 . 
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ELASTIC MODEL 

VISCOELASTIC 

MODEL 

T = Tj / K 


Figure 1. Mechanical block representation of fault. Each block which rests on 
the friction surface is coupled to nearest neighbors and to the driving plate which 
moves to the right with velocity u. (a) elastic coupling clement - perfect spring 
with spring constant K. (b) viscoelastic coupling with spring constants K and K, 
and dashpot viscosity 77 . 
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Figure 3. Event plot for simple elastic force model with 
Kj - 1 x 10 1 7 dyne/cm, Kj = 5 x 10 16 dyne/cm and f* 
randomly distributed in range (2 t 0.03) x 10 20 dyne. 
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Figure 4. Event plot for simple elastic force model with 
Kj ar.d K ; us in Figure 3 but f* randomly distributed in 
range (2 ± 1) x IO 20 dyne. 
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Figure 5. (Continued) 
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Figure 6. Event plot for simple clastic force model with 
Kj and K, as in Figures 3-5 but t* randomly distributed 
in range 0 < f j < 1 x 10 20 dyne. 
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Figure 7. Event plot for viscoelastic force model with 
time dependent friction. The parameters of the sim- 
ulation are = 1 x 10 1 8 , K:j = 3Kj, Kj = 1 x It) 1 7 , 

Kj, =JK jt 1 x 10 :o dyne < f*tt () )<3x 10 20 dyne, ff 
= 0. 98 f; (t 0 ), f. £ (t) = f* (t 0 ) I H 0. 021 log(l + t - t Q )|. 
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Figure 7. (Continu .a) 
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APPENDIX 


SCALING LAWS 


The results of a simulation can be scaled according to the following relations 
derived from Equation 2. 

Given: the results from a simulation performed with independent parameters 

f\ K’, m' 


Define: 


f = Bf' 
K = CK' 
m = Dm' 


Then: 



(A 1) 
(A. 2) 
(A. 3) 

(A. 4) 
(A. 5) 
(A.6) 
(A. 7) 
(A.K) 


where At' is a time interval associated with block motion. We have assumed 


At’ « t;. 
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